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Abstract — The Nimbus I, II, and III satellites provide a new op- 
portunity for climate studies in the 1960s. The rescue of the visible 
and infrared imager data resulted in the utilization of the early 
Nimbus data to determine sea ice extent. A qualitative analysis 
of the early NASA Nimbus missions has revealed Antarctic sea 
ice extents that are significant larger and smaller than the his- 
toric 1979-2012 passive microwave record. The September 1964 
ice mean area is 19.7xl0 6 km 2 ± 0.3xl0 6 km 2 . This is more the 
250,000 km 2 greater than the 19.44xl0 6 km 2 seen in the new 2012 
historic maximum. However, in August 1966 the maximum sea 
ice extent fell to 15.9xl0 6 km 2 ± 0.3xl0 6 km 2 . This is more than 
1.5xl0 6 km 2 below the passive microwave record of 17.5xl0 6 km 2 
set in September of 1986. This variation between 1964 and 1966 
represents a change of maximum sea ice of over 3xl0 6 km 2 in 
just two years. These inter-annual variations while large, are small 
when compared to the Antarctic seasonal cycle. 

Index Terms — Antarctic sea ice, historic data, Nimbus. 


I. INTRODUCTION 

S EA ice extent data from multi-channel passive microwave 
sensors provided one of the longest-running climate 
records from satellites in the form of the sea ice extent index 
[1]. However, while the more than 30-year record is now long 
enough for typical climate-scale change, it is still relatively 
short. Even adding in the Nimbus V single-channel Electrically 
Scanning Microwave Radiometer (ESMR) data extends the 
record by only another 6 years and that data is not as com- 
plete or consistent [2]. Earlier data are only available from 
operational ice charts that are of variable quality and do not 
cover all ice-covered regions. This was particularly true in the 
Antarctic, where there were few ice charts before 1972. This 
paper presents Antarctic sea ice estimates for winter 1 964 and 
1966 derived from early Nimbus satellite data as an initial step 
towards creating Earth Science Data Records of sea ice in the 
1960s. This record has the potential to significantly extend 
the passive microwave record. The Nimbus visible and near 
infrared data, unlike the passive microwave data, sense clouds 
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(very frequent in the polar regions), which may obscure much 
of the surface for days at a time. Satellite visible imagery is 
useful 65° north to 65° south latitude year round. Beyond those 
latitudes, visible light data generally limited to the summer 
periods when there is substantial sunlight. Finally, interpreting 
visible/infrared imagery is more difficult as there are no robust 
automated algorithms like those that have been developed for 
passive microwave imagery [1], [3]. Nonetheless, as discussed 
below, methods have been developed to obtain monthly com- 
posites of sea ice extent from the Nimbus data. 

II. Background 

While early satellite missions during the 1960s collected sig- 
nificant amounts of data, little of it was analyzed at the time. 
The limited use of the 1960s data was generally not due to the 
quality of the raw data, but rather because of the limited pro- 
cessing capability and digital archiving systems available in that 
era. Nimbus I, launched in August 1964, began a sequence of 
experimental satellites built by NASA to develop instruments to 
collect information about the weather (the concept of collecting 
time series satellite climatology was generally not considered at 
the time). Although the instrument and tape recorder only lasted 
three weeks (28 August-22 September) for Nimbus I, it was for- 
tuitous that the measurements occurred in September during the 
likely Arctic ice minimum extent and Antarctic maximum ex- 
tent. The early Nimbus satellites carried the High Resolution In- 
frared Radiometer (HRIR) and the Advanced Vidicon Camera 
System (AVCS) or the Image Dissector Camera System data 
(IDCS). 

The Nimbus HRIR Instrument (Nimbus I, II, and III) used 
a lead Selenide (PbSe) detector to detect infrared radiation to 
a precision of d=l°K in the 3.4- and 4.2-micron bands. The 
Nimbus HRIR data, provides global coverage for the period 
from 1964 to 1972. This newly available data, in conjunc- 
tion with Nimbus IV, V, and 1970s imagery from scanning 
radiometer, infrared and visible sensors on U.S. Department of 
Defense Meteorological Satellite Program (DMSP) platforms, 
could potentially extend the length of the climate record at the 
poles by 50% from the passive microwave 1979-2012 time 
series [4]. Interestingly, the Nimbus I HRIR data has been 
recovered from the original data tapes and likely is the oldest 
satellite-based digital temperature data available anywhere. 
Nimbus II collected data from 15 May 1966 to 18 January 1969. 
Nimbus III collected data from 14 April 1969 to 22 January 
1970. These data acquisition time frames happen to include 
the period when the Arctic sea ice minimums and Antarctic 
maximums occur. Nimbus I did not reach full orbit altitude 
of 1100 km due to a short second-stage burn, resulting in an 
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eccentric 429 km by 937 km orbit [5]. The lower orbit resulted 
in higher spatial resolution data with little overlap between 
images on each orbit. The coverage with Nimbus II in 1966 was 
better because of the higher observational altitude (1100 km) 
and 20% overlap between images, however the data are noisier. 
Nimbus III IDCS images (1100 km) cover a wider extent and 
the images have significantly better dynamic range and overlap. 

The AVCS used video technology to collect snap shots of the 
reflected radiance from the Earth every 91 seconds. The instru- 
ment description can be found in the Nimbus User Guides [6]. 
Navigation of the images is based upon the orbit elements, time 
of measurement and geometry of the AVCS cameras. Although 
most images include a gray scale bar for calibration, variations 
in the camera sensitivity on the satellite and digital-to-analog 
data capture processes make it difficult to construct a “cali- 
brated” radiance from the data. In addition, the data are, at best, 
4-bit resolution so fine details of clouds are not evident. 

III. Methods 

There are three primary steps in obtaining sea ice extent from 
the recovered data: 1) conversion and navigation correction of 
the HRIR digital files recovered by the NASA Goddard Space 
Flight Center (GSFC); 2), digitization, georectification and 
quality control of film images data from AVCS and IDCS held 
by the NO AA National Climatic Data Center; and 3) calculation 
of monthly composite sea ice extents from the recovered data. 

It was critical to recover both the HRIR data as well as the vis- 
ible AVCS/IDCS imagery data since the HRIR provided night- 
time views and AVCS daytime views of the globe. The AVCS 
visible light data is being used to validate the HRIR sea ice ex- 
tents. The value of the HRIR data lies with the fact that during 
the dark polar winter months the HRIR data is the only avail- 
able source. In the case of Nimbus II, the HRIR tape recorder 
lasted nearly three months longer (15 November 1966) than the 
AVCS tape recorder (31 August 1966). 

The GSFC recovered and rescued thousands of Nimbus 
HRIR digital data tapes. The recovery process involved using 
specially-developed tape drives, tape baking, bit detection and 
processing techniques to read the 800 bpi, 7-track tapes [7]. 
The ongoing analysis of HRIR data for sea ice detection is seen 
in Fig. 1. The original Tape Archive Program (TAP format) 
uncorrected (spatially) data from Nimbus I, II and III are 
now available from GSFC. NSIDC will be releasing NetCDF 
corrected versions of this HRIR data in mid-2013. 

The digitization of the AVCS 35 mm film and IDCS 70 
mm film was done at NSIDC. A Kodak HR500 high-speed 
film scanner digitized all the images with 8-bit resolution (the 
original AVCS film was 4-bit) to ensure complete data capture. 
12,000 images from Nimbus I (Fig. 2(a), 2(b)), 120,000 im- 
ages from Nimbus II and 50,000 images from Nimbus 3 were 
scanned. The scanner was calibrated to the grey scale on the 
images. All scanning was performed using consistent fixed pa- 
rameters determined by the calibration to avoid any bias in the 
data capture. A human-augmented custom software application 
documented each image, recording the time and center points. 
This provided navigation accuracy of about ±25 km. A manual 
feature identification technique was developed to identify the 
ice boundary. About 1000 Nimbus I images were reviewed 



Fig. 1. Nimbus II, November 13, 1966, high resolution infrared radiometer 
with possible sea ice edge. 

and an ice edge was identified on about half of those images. 
About 4000 images were reviewed from Nimbus II to create 
four months of ice extents displayed in Table I. 

The Nimbus II data was analyzed by dividing Antarctica into 
2° longitude bins and averaging the observed ice extents for 
each bin. There were many bins with no measured values, es- 
pecially the months of May and June 1966. For short gaps (up 
to 4° longitude) values were interpolated from adjacent values. 
For larger gaps the SSMI climatology was substituted for the 
missing bins. The adjacent interpolation method was used for 
July and August but as there were no large gaps, and substi- 
tuting SSMI values was not required for these months. The per- 
cent complete values for bins are determined before interpola- 
tion or substitution. The ice extent standard deviation is based on 
the average standard deviation of the ice edge latitude in each 
2° longitude bin from all suitable images within the bin. The 
standard deviation of the edge estimates reflect three factors: 
1) spatial variability within the 2° longitude region, 2) temporal 
variability within the month, and 3) uncertainties in the ice edge 
location detected in the images. 

There is considerable correlation between bin measurements 
so it is not reasonable to translate the uncertainty in the bin aver- 
ages directly into the uncertainty of the mean extents. The high 
correlation in bin values leads to significantly fewer than 1 80 in- 
dependent bins. A reasonable assumption is that there are 20 in- 
dependent bin regions with a one-sigma range of area. This un- 
certainty does not include systematic errors, which might arise 
from navigation errors. Qualitatively one can summarize the re- 
sults as saying the May is similar to the SSM/I climatology, June 
and July are larger than SSM/I mean and August is lower than 
SSM/I climatology. 

IV. Error Analysis 

The uncertainties in determining a sea ice extent from this 
early data are primarily tied to the difficulty in distinguishing 
homogeneous clouds from homogeneous ice. Fortunately, ice 
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Fig. 2. a: Sample Nimbus I image September 1964 with a distinct sea ice 
boundary. The red marks show our manual ice edge recordings. These locations 
are used in the average ice extent determination, b: Sample Nimbus I image 
September 1964 with a less distinct sea ice boundary. The red marks show our 
manual ice edge recordings. These locations are used in the average ice extent 
determination. 


edges generally appear different than cloud edges. Additionally, 
clouds tend to move faster than sea ice enabling edge detection 
when utilizing multiple images of the same location. Another 
difficulty in ice edge determination is the lack of independent 
sources of sea ice edge data in the Antarctic during the 1960s. 
The analysis of Arctic sea ice extent agrees very well with the 
historic sea ice charts from the 1960s (where the data overlaps) 
giving further credibility to the project results. 

In an effort to create an automated ice detection method, a 
determination of ice extent was done by remapping and com- 
positing the minimum brightness for all the images to derive a 
clear-sky radiance. This gives a qualitative impression of the lo- 
cation of the ice edge and provides an independent validation of 


the ice extent [8], but because of variations in the calibration, the 
project did not succeed in using the composited image for an au- 
tomatic ice edge analysis. For the Nimbus II data a similar veri- 
fication technique was developed. An estimated ice edge is cal- 
culated from a running 1 1 -day minimum brightness composite 
for the images with an obvious sea ice edge during August 1966 
(Fig. 3). Visually, it shows a solid match between the average 
manual ice edge and the composite minimum brightness image. 
There are some pixels with anomalous brightness (noise) in the 
image, which is why the method is not suitable for an automatic 
ice detection scheme, but it is useful for a qualitative evalua- 
tion of the manual analysis. It may be possible to minimize the 
noise by modifying ice edge techniques developed for synthetic 
aperture radar data with high speckle noise [9], allowing a di- 
rect estimate from the minimum brightness composite, but that 
is beyond the scope of this study. 

V. Antarctic Sea Ice Extent Estimates, 1964 vs. 1966 

The focus of this paper is on the comparison of anomalous 
maximum sea ice extents from September 1964 (Nimbus I) and 
August 1966 (Nimbus II) using AVCS visible imagery. Early 
analysis of the 1969 AVCS data indicates a return to less anoma- 
lous conditions. 

A. Antarctic Extent 1964 

Antarctic sea ice extent for 1964 was calculated from the 
Nimbus I imagery, as described in [10]. The 1964 satellite data 
provide fairly complete coverage over the Antarctic, allowing 
the ice edge to be clearly seen in most areas in at least one image 
over the three weeks. This provides a reasonably robust estimate 
of sea ice extent for the period, though some uncertainty results 
from the variation of the ice edge location during the time pe- 
riod. The Nimbus I data had at most eight images of any area 
in the three weeks of data collection, so it was difficult to dis- 
tinguish homogeneous clouds and homogeneous ice. As illus- 
trated in Fig. 4, the manual sea ice analysis focused on the ice 
margin. The analyzed ice edge locations were used in the av- 
erage ice extent determination, making the assumption that lo- 
cations continent- ward of the edge were ice-covered (i.e., >15% 
concentration sea ice). Comparing modern SSM/I and AVHRR 
visible light records could further test this assumption, however 
this was beyond the project scope. 

Fig. 5 shows the mean and minimum and maximum ice edges 
for September 1964, from an average in 2° longitude bins. That 
is, all the individual latitudes falling in a longitude bin were 
averaged with the range recorded. Missing bins were filled with 
an interpolation from nearby good measurements to make the 
full area average up to 4° . For larger gaps, as found in May and 
June 1966, SSM/I climatology values were substituted for the 
missing values. For comparison the SSM/I 1979-2000 median 
climatology is plotted. Qualitatively, the Nimbus ice edge is 
often north of the passive microwave results [4] suggesting 
more ice than the climatology. The ice area inside the mean 
edge is 19.7 x 10 6 km 2 ± 0.3 x 10 6 km 2 . The uncertainty is 
based on the difference between the minimum and maximum 
area observed. Areas of the small polynyas were not subtracted, 
to better match the lower resolution of the passive microwave 
record, which can only resolve very large open water areas 
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TABLE I 

Antarctic 1966 Sea Ice Extent Based Manual Analysis of Nimbus AVCS Visible Images. Estimates From Images Were Composited Into a Single 
Extent Estimate Within a 2° Longitude Bin (180 Bins Total). The Bins Were Classified as Having “High Quality” Images (i.e., Images Where 
the Ice Edge was Clearly Apparent) and “Medium Quality” Images (i.e., Images Where There Were Ambiguities About the Ice Edge Locations). 
Bins That Did Not Have Images or Only had “Low Quality” Images Were Filled With Interpolated Edge Estimates From Surrounding Bins 


Month 

Nimbus Antarctic Ice 
extent based only 
on high quality bins 
and interpolation (± 
1 st. dev. range 
within bins) 

Nimbus Antarctic 
Ice extent based on 
high and medium 
quality bins and 
interpolation (± 1 
st, dev. range 
within bins) 

% of bins 
with high 
quality 
edge 

estimates 

Ice edge 
standard 
deviation 
for high 
quality 
bins 

1979-2000 SSM/I 
climatology ice 
extent and one st, 
dev. range over the 
climatology period 

May 

N/A 

10.7 ±0.2 x 10 6 km 2 

N/A 

N/A 

10.3 +0.56 x 10 s km 2 

June 

14.3 +0.4 xlO 6 km 2 

14.8 +0.3 x 10 s km 2 

75 

0,81 

13.3 +0.40 xlO 6 km 2 

July 

16.4 ±0.5 x 10 6 km 2 

16.4 ±0.4 x 10 s km 2 

98 

1-05 

15.8 ±0.30 xlO 6 km 2 

August 

15.9 +0.4 xlO 6 km 2 

16.1 +0.4 x 10 6 km 2 

99 

0,76 

17,5 +0.22 xlO 6 km 2 



Fig. 3. Minimum AVCS brightness for the last week of August 1966 for images 
with obvious sea ice edges. Superimposed is the sea ice edge (black line) from 
the manual analysis. 

within the ice edge. Notice that the 30°E and 90°E scallops 
appear in both data sets. There are no corroborating ship or 
aircraft measurements in 1 964 so this is a stand-alone measure- 
ment that cannot be directly validated against comparison data. 
However, the estimate is consistent with Predoehl 1966 [11], 
which was estimated from the same Nimbus I raw data. 

To put this ice extent number in context, the time variation of 
the passive microwave observations [1] from 1979 to 2010 are 
shown in Fig. 6. The range of the passive microwave in each 
year is the mean ± 1 standard deviation of the daily area values 
for each September. September 1964 Nimbus I sea ice estimate 
is the largest Antarctic sea ice extent ever observed by satel- 
lites [12]. September 1964 Arctic extent was also estimated by 
Meier et al. [10]. The estimate is near the 1979-2000 clima- 
tological average and analyzed ice edge locations are consis- 
tent with other Arctic estimates in 1964 based on operational 
ice charts. 

B. Antarctic Extent 1966 

Coverage of Nimbus II AVCS images for the Antarctic for 
1966 was fairly complete but the film record ends on August 
31 so a direct comparison with the 1964 estimate from Nimbus 


■ SSMI Climatology Mean Nimbus Range at each la ngatude 



Fig. 4. Mean Antarctic ice edges Nimbus I with the range of minimum and 
maximum observed ice extent. The grey dashed line is the SSM/I climatology 
from 1978 to 2000. 

I is not possible. The manual analysis procedure for AVCS im- 
ages was improved by establishing a quality indicator (2 = best, 
1 = usable but uncertain) for the ice edges. Fig. 6 shows the lo- 
cations of ice edges for the whole month of August 1966 from 
the best and useable quality locations. This August 1966 data 
is then compared to the passive microwave 1979-2000 clima- 
tological ice extent for August. The results for May 1966 had 
higher uncertainty because images often lacked clean ice edges. 
For both May and June there were longitudes with no ice edge 
detection, leading to further uncertainty. The ice extent esti- 
mates for June, July and August, based on images with clean ice 
edges, are 14.6 x 10 6 km 2 , 16.4 x 10 6 km 2 and 15.9 x 10 6 km 2 , 
respectively. The 1966 June and July ice extents are just slightly 
below the minimal historic values for the sea ice extent relative 
to the passive microwave record [4], [12], [13]. However, the 
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Fig. 5. Comparison of SSM/I September Antarctic ice extent with Nimbus I 
1964 result. The ranges for the SSM/I are the 1 standard deviation of the 3-day 
means. The extremes for the 1964 data are the averages of all the maximums 
or minimums in each longitude bin. The box range represents d= one standard 
deviation each bin average. 


SSMI -"ICE EDGE AVERAGE EDGE 



Fig. 6. Map of Antarctic ice edges for August 1966 with SSM/I climatology in 
green. Each red dot represents one ice edge determination. 

August 1966 estimate fell well below the passive microwave 
record’s minimum range. The Nimbus II ice extent of 15.9 mil- 
lion km 2 is 2.2 million km 2 lower than the passive microwave 
1979-2000 average of 18.1 million km 2 , and is 1.7 million km 2 
lower than any August extent from 1979 to 2012 (Fig. 6). It is 
also lower than any individual daily values in the passive mi- 
crowave time series for August. The anomalously low extent 
is due to the Nimbus ice edge being observed further poleward 
than in the passive microwave data, especially near 170°E. 



Fig. 7. Comparison of SSM/I August Antarctic ice extent with Nimbus II 1 966 
result. The range shown for the SSM/I data (1979-2013) is the standard devia- 
tion of the daily areas. 

The fact that our August 1966 is such an outlier indicates po- 
tential issues in the Nimbus August data. In addition, the drop in 
extent from July to August is not seen in the passive microwave 
record and is not physically likely. This issue may be related to 
the data quality, the effect of the timing of the imagery on the 
monthly composite estimate, limitations in the manual analysis, 
or other factors. AVCS data is not available for September, the 
month when the Antarctic typically reaches its maximum. On 
the other hand, in the 1972 analysis of the 1966 Nimbus II data, 
Sissala [14] reported an even lower August ice extent of 15.17 
million km 2 . 

One potential explanation comes from the intra-month vari- 
ability of the ice cover during August. Notice in Fig. 7 that there 
is significant variation in passive microwave ice extent in Au- 
gust, much more than in September. This is because the ice 
extent is still increasing in August as freezing occurs. So tem- 
poral sampling gaps or biases (i.e., sampling more early in the 
month vs. late in the month) can potentially bias a monthly com- 
posite estimate. Near the time of maximum extent in September 
both freezing and melting are occurring so there is less large- 
scale daily change and the potential for a sampling bias will be 
smaller. Qualitatively, we do see the high daily variability the 
Nimbus II data for August 1966. There is not significant vari- 
ability in the Nimbus I 1964 data in September, due in part, to 
the shorter 3 -week record. 

Unfortunately, we have been able to uncover any corrobo- 
rating information on Antarctic sea ice in 1964 or 1966, in- 
cluding any data or imagery from US Navy aircraft over-flights 
[15]. There are a few other potential sources of satellite data 
from this time. The Soviet satellite Kosmos 44, 45 and 46 mis- 
sions flew in August and September of 1964 [16], however 
the continuing existence of this data is in doubt. There was 
considerable DMSP data captured in the 1960s. Unfortunately, 
the Air Force indicated that all the DMSP imagery from the 
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Fig. 8. August 1 964, Nimbus I Antarctic AVCS image of the same general area 
as the Wendell Polynya seen the 1974-1976 Nimbus V ESMR data. Note the 
abundance leads and fractures. 


1960s was destroyed and reclaimed for its silver content after it 
was analyzed (2010, personal communication, Dennis Hobson, 
DAFC Deputy Chief, Operational Capability Requirements Di- 
vision HQ AFWA/A5R). Two 1964 Argon K5 satellites col- 
lected data from June-July and August-September 1 964 respec- 
tively, but only a portion of the data is currently available [17] 
and will require further investigation. A 1 963 Argon K5 dataset 
was recently reprocessed to produce an image of the Antarc- 
tica continental edge. This data may have potential for sea ice 
extent information as well [18]. The Wendell polynya seen in 
the 1974-1976 ESMR data did not appear in the Nimbus I or 
II data, however the 1964 data reveals numerous fractures and 
leads seen the same area as this polynya (Fig. 8). 

VI. Summary 

This study presents a method to calculate a composite 
monthly sea ice extent from the AVCS visible sensor on 
Nimbus satellites that were in orbit starting in 1964. Esti- 
mates from Nimbus I and Nimbus II for September 1964 and 
May- August 1996 respectively found anomalous Antarctic sea 
ice extents in both 1964 and 1966. The extent of sea ice around 
Antarctica was larger in September 1964 than in any September 
measured since (1973-1976 and 1979-2012) [1], [2], [11]. 
Even including the error bar, the 1964 results are outside the 
range of recent data, though 2012, a record high in the passive 
microwave record is on the outer range of the Nimbus estimate 
error bars. In 1966, extent was within 1973-2012 levels during 
May-June, but the August 1966 Antarctic extent was lower 
than any other August in the satellite record [4], [12], [13]. 
This potentially indicates an issue with August 1966 imagery. 
However, qualitative evaluation of the imagery and previous 
analysis of the imagery [14] appear to verify our estimate. Even 
accounting for the seasonality difference of 3.8 x 10 6 km 2 
between August and September, the sea ice change between 
winter 1964 and 1966 is unprecedented in the Antarctic climate 
record. To put this in perspective, the unusual decrease in sea 
ice extent between 1973 and 1976 was only 2.17 x 10 6 km 2 


[12]. Future analysis of the HRIR data may yield an estimate 
for September-November 1966 for improved comparisons. 

In addition to the Nimbus data, we are attempting to recover 
the polar orbiting TIROS 9 and 10 visible data from film from 
1965 and 1966 to fill in some temporal gaps. Scanning of the 
data is straight forward, but the documentation and metadata is 
very limited so navigation of the images will be difficult. These 
film images do not contain any time stamps; however the time 
of the start of the orbit is generally available so interpolation 
with selected feature identification may enable navigation that 
is accurate enough for the data to be useable. 

In the future, an analysis using the same methods developed 
for this study on post- 1973 DMSP visible data and comparing it 
with Nimbus 5 ESMR data from the same temporal range will 
further test the validity of this methodology. As seen from the 
sea ice examples, preservation of the early Nimbus data can 
be a valuable climate data record for observing important cli- 
mate change indicators. In addition to extending the length of 
the 1979-2011 passive microwave sea ice records by as much 
as 50% (16 years), this data could be used for characterizing 
the Earth’s climate and model validation for other parts of the 
globe including: sea surface temperatures, cloud fraction, hur- 
ricane, and potentially some land surface characteristics. These 
data will be available in NetCFD format from NSIDC in the 
summer of 2013. 
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